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Abstract. Strong and fragile glass relaxation behaviours are obtained simply 
changing the constraints of the kinetically constrained Ising chain from symmetric 
to purely asymmetric. We study the out-of-equilibrium dynamics of those two models 
focusing on the Kovacs effect and the fluctuation-dissipation relations. The Kovacs or 
memory effect, commonly observed in structural glasses, is present for both constraints 
but enhanced with the asymmetric ones. Most surprisingly, the related fluctuation- 
dissipation (FD) relations satisfy the FD theorem in both cases. This result strongly 
differs from the simple quenching procedure where the asymmetric model presents 
strong deviations from the FD theorem. 
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1. Introduction 

The Kovacs or memory effect has been observed by Kovacs himself in the 1960s on 
structural glassy systems pQ. It is a surprising memory effect of the energy (or volume) 
of the system following a particular quenching procedure. This quenching procedure 
consists in suddenly cooling a glassy system from a high temperature (infinite one in 
our case) to a very low intermediate temperature Tj and letting the system relax. When 
the waiting time t w necessary for the energy of the system to reach the equilibrium value 
for a final temperature Tf is attained (e(t w ) = e eq (Tf)), the temperature of the system 
is set to this value Tf > T. The result found by Kovacs is that, even though the system 
is at the equilibrium energy (the volume in his case) corresponding to the temperature 
imposed, the system and its energy are still evolving. The system keeps memory of its 
history and of the fact that equilibrium is not effective. After a rapid increase of the 
energy, the system reaches the equilibrium and the energy decreases and levels off that 
equilibrium leading to a hump in the energy as function of time. 

In this paper, we are interested in a comparison of the Kovacs effect for two 
simple (even simplistic) models with respectively strong and fragile glass behaviours. 
We consider the symmetric and purely asymmetric kinetically constrained Ising chain 
(KCIC) models j2JE]- This simple change of the kinetic constraints allows one to switch 
respectively from strong to fragile glass behaviour. The underlying equilibrium being 
the same for both models, a direct comparison of the dynamical effects is possible. 

After a short presentation of the models in section 2, the Kovacs effect is discussed 
in section 3. The effect is observed in both models considered but is enhanced in 
the asymmetric case. This section also contains simple rescaling arguments to explain 
this effect and a comparison with recent works on the Kovacs effect [U El EJ E|- 
The fluctuation-dissipation (FD) relations are studied during the Kovacs quenching 
procedure and presented in section 4. The FD relations satisfy the FD theorem in 
both cases. These results are in strong contradiction with those obtained using a 
simple quenching procedure from high temperature to a low final temperature after a 
waiting time t w . This last quenching procedure leads to FD relations satisfying the FD 
theorem for the symmetric KCIC model and for waiting times well below the relaxation 
time. In contrast, strong deviations from the FD theorem have been observed for the 
asymmetric KCIC model for waiting times smaller than the equilibration time. We give 
some conclusions in section 5. 

2. Presentation of the models 

In this paper, we are interested in the possible difference concerning the Kovacs effect 
due to strong or fragile glass behaviour. We thus consider the KCIC model [21 El for 
which constraints may be chosen to model strong (for symmetric ones) and fragile (for 
a purely asymmetric chain) glass relaxations. 

Let us consider a chain of iV Ising spins (o"j = 0,1 with i = 1,---,N) without 
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interactions where spins <7j = 1 are considered as defects. The corresponding 
Hamiltonian is thus trivial (H = ^\ af) as well as the equilibrium thermodynamic 
properties. The equilibrium energy at temperature T or inverse temperature (3 is given 
by e eq (T = 1//3) = 1/ (1+e 13 ) which is also the concentration of defects or the probability 
to have a defect at site i. It is possible to determine exactly the probability for a defect 
to have its next defect (on the left) at a distance d 



The first term on the right hand side of the equation corresponds to the probability to 
have a defect whereas the second term (with a power d — 1) is the probability to have no 
defects in the intermediate d — 1 sites. The non-interacting spins render at equilibrium 
the probabilities at each sites independent of each other and leads to this simple product 



All these equilibrium properties are independent of any dynamics considered. 
However, the introduction of kinetic constraints allows one to obtain a slowing down 
of the dynamics characteristic of glassy systems before the equilibrium properties are 
reached. The probability for a spin to flip is constrained in the following way: in the 
symmetric case, a spin is able to flip as soon as a neighbour (left or right) is a defect 
whereas, in the asymmetric case, the defect has to be on the left. Such spins are 
also called spin facilitated and their probability transitions are given by the following 
equation: 



The first term on the right hand side corresponds to the usual Metropolis probability and 
allows one to satisfy the detailed balance. The second term corresponds to the general 
kinetic constraints with a probability b to flip the spin if there is a left neighbour and 
1 — b for a right neighbour. The symmetric model (b = 1/2) and the purely asymmetric 
one (6 = or 1) correspond to particular values of this parameter b. 

These models have been extensively studied (for more information and references 
on kinetically constrained models see the recent review by Ritort and Sollich 8\). With 
the symmetric constraints, the dynamical behaviour is reminiscent of a strong glass 
with a relaxation time following an Arrhenius law. At sufficiently low temperature, 
the defects are mainly isolated and may be considered as simple particles diffusing 
with a temperature-dependent rate of diffusion T ~ exp(— 1/T). The energy (or 
concentration of particles) evolves through creation and annihilation processes. Similar 
react ion- diffusion models have been introduced for a long time to study domain growth, 
coarsening and aging [HI E3 E] • Within the asymmetric constraints, the energy barriers 
involved in the motion of defects are increasing logarithmically with the distance from 
the next defect [12]. As a consequence, the relaxation time follows the Bassler law [13J: 
treiax ~ exp(l/T 2 In 2) . Whereas the Arrhenius behaviour is associated to a strong glass 
behaviour following the Angell's classification Jl] , the super- Arrhenius behaviour of the 
asymmetric model is associated to a fragile glass. Intermediate constraints allow the 




(1) 



in P eq (d,T). 
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Figure 1. Left : Kovacs effect for a final temperature Tf = 1/3 and different 
intermediate temperatures Tj within the symmetric KC1C model: the energy as 
function of time elapsed since the waiting time is plotted. Right: the rescaled energy 
Seif for short times (see text for definition) is plotted for two final temperatures 
Tf = 1/2 (full lines) and 1/3 (dashed lines) and for different intermediate temperatures. 



system to continuously crossover from fragile to strong glass behaviour P2] but will not 
be considered in this study. 

3. Kovacs effect 

As already mentioned in the introduction, the Kovacs effect is observed following 
a particular quenching procedure. At time t — 0, a system, equilibrated at high 
temperature (T = oo in our case), is suddenly quenched to an intermediate low 
temperature T«. The system starts to relax and the energy decreases until it reaches the 
equilibrium energy corresponding to a final temperature Tf > Ti after a waiting time 
t w (Ti,Tf) defined by the following equation: 

e (*«) = ^$>^») = ee ^ T /) = (! + e/3/ ) _1 (3) 

i 

where e(t) is the energy of the system at time t and f3 — 1/T is the inverse temperature. 
The temperature of the system is set to Tf at this waiting time t w . 

If the system was characterized only by the thermodynamical parameters (energy, 
volume and temperature), we would expect the energy of the system to stay constant and 
equal to e eq (Tf) after the waiting time t w . However, even though its energy corresponds 
to that equilibrium one at the imposed temperature, the particular configuration of the 
system at t w is still far from an equilibrium configuration at Tf. As a consequence, the 
energy is still evolving after t w . 
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3.1. The symmetric KCIC model 

The figure H illustrates the Kovacs effect for the symmetric KCIC model ^Hj . As can 
be seen on the left panel, the Kovacs effect is an increasing function of the difference 
between the intermediate and final temperatures. A two timescale behaviour is also 
clearly observed. The energy increases on a timescale t ~ 0(1) to reach a plateau 
followed by a further increase up to a maximum. The energy finally decreases to reach 
the equilibrium energy at temperature Tf. Note also that the maximum in energy 
decreases and shifts to slightly higher times for higher intermediate temperatures Tj. 
This result was already observed for the volume by Kovacs [T| and for the energy in 
recent simulations |U E] . 

The maximum of the energy is obtained for a timescale t m ~ 0(exp(3/Tj)). 
This timescale is of the same order as the equilibrium timescale t eq (time to reach 
the equilibrium after a rapid quench from a high temperature to a low temperature 
Tf) jHlEj- This exp(3/T/) dependence of t eq is explained considering that the energy 
e(t) after such a quench decays like the annihilation process A + A — > A in one 
dimension with a diffusion rate T: e(t) ~ (rt)~ 1//2 . The energy thus reaches that 
equilibrium (e eq (Tf) ~ e _1 / T/ ) after an equilibrium time t eq ~ e~ 2 ~ e 3 ^ Tf due to 
the temperature-dependent diffusion rate (r ~ e -1//Tf ). With a similar argument it is 
possible to estimate the waiting time t w (Ti, Tf) ~ r _1 (Tj) e~ 2 (T/) ~ e^ Ti e 2 ^ Tf where the 
equilibrium energy to reach is e eq (Tf) ~ e~ l l Tf but the diffusion rate is T(Tj) ~ e~ x l Ti 
since the temperature of the system before t w is set to Tj. We recover the equilibrium 
time when Tj = Tf. The estimate of the waiting time agrees qualitatively with the 
values obtained from numerical simulations for the different couples of intermediate and 
final temperatures considered. 

Let us now analyse the short timescales behaviour. The fast increase of the energy 
is related to the probability pi to have neighbouring defects. Such a probability is 
given at t w by pi = e(t w ) e eq (Ti) where e(t w ) = e eq (Tf) is the probability to have a 
defect at t w and e e(J (Tj) is the probability for its neighbour to also be a defect. This 
last term corresponds to the equilibrium probability P eq (d = l,Tf) due to the fact 
that the timescale to equilibrate the concentration of neighbouring defects (or non- 
constrained defects) is t ~ 0(1) <C t w and has thus already equilibrated at t w . When 
the temperature is increased from the low temperature Tj to a higher one Tf, on a 
similar timescale t — t w ~ 0(1) independent of the temperature, the probability to 
have neighbouring defects reaches its equilibrium value Pf = e 2 eq (Tf). Furthermore, the 
timescale for the motion of defects is r _1 (Tj) ~ exp(l/Tj) ^> 1. As a consequence, the 
change of probability from pi to pf is mainly obtained from the creation of defects and 
the energy as function of time may be expressed as 

e(t - t w ) = e eq + (pf - Pi)5e if (t - t w ) (4) 

where we expect 5eif(t — t w ) to be independent of both intermediate and final 
temperatures for t — t w ~ 0(1) and to start to differ at <5ej/ ~ 1 for different Tf. This 
prediction is verified on the right panel of figure ^ where 5 ej / is plotted for different 
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Figure 2. Left : Rescaled energy Ae,/ as function of the rescaled time (t—t w )/t m (see 
text for definition) and the distribution of distances between defects at the different 
waiting times t w (Ti,Tf). For both figures, the model considered is the symmetric 
KCIC and the final temperature is Tf — 1/3. 



couples (Tj, Tf) on short times t — t w . 

On longer timescales, the activated regime is involved. The defects may then 
be considered as simple particles moving with a temperature-dependent diffusion rate 
Tf ~ exp(— 1/Ty). The number of particles (or energy) evolves through creation of 
neighbouring defects and annihilation when two defects collide (A + A <-> A). As we 
have seen with the short timescales behaviour, the energy increases after the waiting 
time t w . This increase is not restricted to the short timescales but continues on longer 
timescales and is due to the fact that the distribution of distances between defects P(d) 
at the different waiting times t w (Ti,Tf) is far from the equilibrium one P eq (d,Tf) (see 
Eq. © for definition) . This difference is evident on the right panel of figure El The 
distribution P(d) converges to that equilibrium when the intermediate temperature 
reaches the final one as it should. One important point is that the whole range of the 
distribution is affected. Defects with neighbouring defects at all distances relax with 
the same typical timescale due to the single energy barrier involved (the motion of a 
defect from one site to the next one is independent of the distance from its neighbouring 
defect). 

The interpretation with diffusive particles and with annihilation-creation processes 
is similar to coarsening models. The only differences with the ID Ising model concern the 
diffusive rate which is temperature-dependent in our case and the creation-annihilation 
processes which are of the type A + A <-+ in the Ising model. This analogy to domain 
growth models allows us to use the same rescaling for long timescales of the Kovacs effect. 
In recent works jHEHH], it has been shown that, for domain growth models, the energy 
shift from the equilibrium energy may be rescaled for all intermediate temperatures T 
in the following way: 




(5) 
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Figure 3. Comparison of the Kovacs effect for the symmetric (dashed lines) and the 
asymmetric (full lines) KCIC models. The energy as function of time elapsed since 
the waiting time is plotted for a final temperature Tf — 1/2 and different intermediate 
temperatures Ti = 1/4, 1/3 and 1/2.5. 



where t m is the time for which the shift in energy from the equilibrium energy already 
diminished by a factor two from its maximum value 

e(t m ) = e eq (T f ) + e m /2. (6) 

As can be seen on the left panel of figureEl this rescaling is correct for the long timescales. 
For short timescales, a small deviation from this rescaling is observed. 

3.2. The asymmetric KCIC model and comparison 

In figure El we compare the Kovacs effects for the symmetric and asymmetric KCIC 
models. The energy as function of the time t — t w since the final temperature Tf = 1/2 
was set to the system is plotted for different intermediate temperatures Tj = 1/4, 1/3 and 
1/2.5. As can be seen, the short timescales behaviour (non-activated) is similar in both 
models. The same rescaling introduced for the symmetric KCIC model works perfectly 
on short timescales for the asymmetric KCIC model (see left panel of figure HJ. The 
argument developed to obtain the rescaling is independent of the dynamical constraints 
as soon as the motion of particles occurs on timescales larger than the spontaneous 
creation of neighbouring defects. 

On longer timescales, when activation plays an important role, the Kovacs effect 
starts to differ between the symmetric and asymmetric KCIC models. The main result 
is a higher maximum of the energy for the asymmetric case. The timescales involved 
are also not surprisingly different due to equilibration times which strongly differ from 
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Figure 4. Kovacs effect for the final temperature Tf = 1/2 and different intermediate 
temperatures T, with the asymmetric KCIC model. On the left figure are plotted the 
rescaled energy Seif for short timescales as function of time since the waiting time and, 
on the right figure, the rescaled energy for long timescales as function of the rescaled 
time. 



strong to fragile glass behaviour. This can also be traced out on the values of the 
waiting times t w (Ti, Tf) larger for the asymmetric KCIC model (fragile glass behaviour) 
than for the symmetric KCIC model (strong glass behaviour). From the analysis of 
Sollich and Evans [Ej, it is possible to have an estimate of the waiting time t w (Ti,Tf) 
for the asymmetric KCIC model. The same creation-annihilation processes occur as 
for the symmetric KCIC model but the motion of particles may not be considered as 
simple diffusion due to the increasing energy barriers with the distance between defects. 
An anomalous coarsening occurs where the energy decreases like e(t) ~ ^-Tiin2 with a 
temperature-dependent exponent in contrast to normal coarsening. The waiting time 
is obtained when the energy reaches the equilibrium energy e eq (T f ) ~ e- 1 '^ leading 
to t w (Ti,Tf) ~ exp(l /TiTf In 2). Note that once again we recover the equilibrium 
timescale t eq (Tf) ~ exp(l/T|ln2) when Tj = Tf. The estimate for the waiting time 
agrees qualitatively with the values obtained from the numerical simulations for different 
couples of intermediate and final temperatures. 

We use the rescaling for the long timescales of coarsening models (JHJ (see right 
panel of figure 0} as for the symmetric KCIC model. We observe some deviations from 
a perfect rescaling even for the long timescales. This difference is due to the anomalous 
coarsening behaviour of the asymmetric KCIC model. The failure of the rescaling on 
short timescales is also more evident. 

4. Fluctuation-dissipation relations 

The violation of the fluctuation-dissipation (FD) theorem is usually present in glassy 
systems due to the out-of-equilibrium dynamics (see for example the recent review 
by Crisanti and Ritort ^H] and the references therein). FD relations have already 
been considered for the symmetric and asymmetric KCIC models in jTH] with the 
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Figure 5. Fluctuation-dissipation relations for the symmetric (left) and asymmetric 
(right) KCIC models with a usual quenching procedure (see text for explanation). A set 
of final temperatures (X/ = 1//?/ = 1/3 and 1/6) and waiting times (t w = 10, 100, 1000 
and 10000) have been studied and all couples of temperatures-waiting times (dotted 
lines) satisfy the FD theorem (full line) in the symmetric case. A set of final 
temperatures (Tf = l/(3f = 3 and 5) and waiting times (t w = 10 and 1000) have 
been studied and strong deviation from the FD theorem (full line) have been observed 
for all couples of temperatures-waiting times in the asymmetric case. 



usual quenching procedure. A more careful analysis of the symmetric case leads to 
the conclusion that even far from equilibrium (at least after a short transient), the 
FD relations satisfy the FD theorem for all temperatures and waiting times [20] ■ This 
result was explained by the fact that defects are mainly isolated at low temperatures 
after a short transient and behave as in local equilibrium (with a temperature-dependent 
diffusive rate T ~ exp(— 1/T)) even if the number of defects is larger than in equilibrium. 
This argument is no longer valid for the asymmetric case where the motion of a defect 
is dependent on the position of its left neighbour defect and is thus sensitive to the total 
number of defects or at least to the distribution of distances between defects which 
is out-of-equilibrium. The differences between the FD relations for both models are 
evident in figure El following the usual quenching procedure. We consider such FD 
relations after the Kovacs quenching procedure to check if this discrepancy between the 
two models is also present. 

Let us first explain how those FD relations are determined. The integrated response 
is determined following the now standard procedure introduced by Barrat |21j . A 
pertubation 5H(t) = —h(t) EiOi is added to the Hamiltonian where hit) = h Q(t—t w ) 
is a field of strength h introduced after the waiting time t w (6(x) is the Heaviside 
function). Random fields are set to each site through the random variables £j = ±1 
with identical probabilities. The integrated response is then 

X(t,t w ) = -^^2£i(ai(t)) h (7) 

i 

with the overline standing for an average over the random variables £j and the 
brackets for the dynamical average in the presence of the perturbation. Note that 




Figure 6. Fluctuation-dissipation relation for the symmetric (left) and asymmetric 
(right) KCIC models with a Kovacs quenching procedure (see text for explanation). A 
set of intermediate and final temperatures (Tj = l//3j, Tf — 1//3/) have been studied 
and all couples of temperatures (dotted lines) satisfy the FD theorem (full line). Final 
and intermediate inverse temperatures considered for the symmetric case: /?/ = 1,2 
and 3 and 0i = 3,4 and 5 and for the asymmetric case: f3f — 1,2 and 3 and /?, = 3 
and 4. 



in these kinetically constrained models different dynamics may be considered after the 
introduction of the perturbation. We considered the modified Metropolis one discussed 
in [22] which gives similar results to the Metropolis one in this case 

The integrated response has to be compared to the corresponding correlations to 
check the validity of the FD theorem. In the case where the energy is still evolving the 
correct correlations to consider are the connected ones 

C(t, t w ) = - J^Oift)^)) - e(t)e(t w ) (8) 

i 

where e(t) is the energy at time t and the brackets stand for the dynamical average 
without the perturbation. In equilibrium, the integrated response and correlations only 
depend on the difference of times and satisfy: Tx(t — t w ) = C(0) — C(t — t w ). It has 
been shown in [20J that a correct generalization of this expression for systems out-of- 
equilibrium is 

T X (t,t w ) =f(Y(t,t w )) (9) 

with Y(t,t w ) = C{t,t) — C(t,t w ). If the FD theorem is satisfied the function f{x) = x. 
Thus, a parametric plot of the integrated response with respect to the function Y allows 
one to check the validity of the FD theorem. Note that in these kinetically constrained 
systems, C(t,t) = e(t) — e 2 {t) depends on time. This explains the parametric plot with 
respect to Y(t,t w ) instead of the two-time correlation C(t,t w ). A slope X different 
than 1 in this parametric plot would suggest the existence of an effective temperature 
T eff = T/X. 

In figure El we have plotted the FD relations for both symmetric (left) and 
asymmetric (right) KCIC models following the usual quenching procedure. At t — 0, 




Figure 7. Responses (full lines), energies (dotted lines) and correlations (dashed lines) 
for different intermediate temperatures Ti (from top to bottom Ti = 1/5, 1/4 and 1/3) 
for the symmetric KCIC model with Tf = 1/3. 



the system is quenched from an infinite temperature to a low temperature Tf. The 
perturbation is introduced at a waiting time t w after this quench. On the left panel, the 
FD relations for the symmetric KCIC model satisfy the FD theorem for a large set of 
temperatures and waiting times (as soon as t w > 1). The validity of the FD theorem is 
far from obvious if we keep in mind that both the integrated response and Y(t,t w ) are 
non-monotonic. However, they compensate each other to follow the FD theorem curve 
at all times. The non-monotonic behaviour of the response is due to the decrease of 
the energy on long timescales. When rescaled by the energy, the response is continously 
increasing. 

For the asymmetric KCIC model, the FD relations show strong deviation from the 
FD theorem. The non-monotonic integrated response cannot be accounted for by the 
difference of correlations Y(t,t w ) as it is for the symmetric KCIC model. It leads to 
different humps in the FD relations corresponding to the different plateaus observed in 
the decay of the energy. The increasing energy barriers (AE = 1,2,...) lead to those 
different plateaus in the energy decay after a quench in temperature [T2"l ITS] . The FD 
theorem is recovered for sufficiently large waiting times (for t w > 10 4 when Tf = 1/3 for 
example). 

Let us now consider the FD relations for the symmetric (left) and asymmetric (right) 
KCIC models after the Kovacs quenching procedure (see figure EJ)- As can be seen in 
both cases the FD theorem is verified. This is in strong contrast with the evolution 
of the energy which shows a deviation from equilibrium. This deviation is sufficiently 
important to exclude the fact that it does not show up on the FD relations. In figure [7| 
we can observe that the responses for different intermediate temperatures as well as 
the energies and the correlations are different. However, the FD theorem is satisfied. 
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Note for example that when T$ is highly slower than Tf, the response shows a hump 
which may be related directly to the fact that the energy is not constant. An increasing 
energy, corresponding to a larger number of defects, tends to increase the response. The 
further decrease of the response is then related to the same decrease in the energy. The 
consequence on the FD relations is that the parametric curve increases from the origin 
to a maximum before a small decrease occurs but always following the line f(x) = x 
corresponding to the FD theorem curve. 

The fact that the FD relations follow the FD theorem curve for the symmetric KCIC 
model is not surprising since it is already the case after a usual quenching procedure. 
A simple explanation suggesting that defects are in local equilibrium has been given. 
The surprise comes from the asymmetric KCIC model where the FD relations differ 
strongly from the FD theorem curve for most waiting times with the usual quenching 
procedure. With the Kovacs one, the FD relations follow the FD theorem curve. A 
possible argument would be the following. The waiting times involved in the Kovacs 
quenching procedure are sufficiently large that the out-of-equilibrium configuration may 
be considered as equivalent to a simple fluctuation in the system at equilibrium. In 
other words the linear regime is attained and the FD theorem is satisfied. However, this 
argument would not explain the shift in energy observed. 

5. Conclusion 

In this paper, we have studied the Kovacs effect for the symmetric and asymmetric KCIC 
models with respectively strong and fragile glass behaviours. The effect is stronger 
in the asymmetric case. The symmetric KCIC model satisfies the same rescaling as 
domain growth models for long timescales. Some small deviations are present for the 
asymmetric KCIC model and may be explained due to the anomalous coarsening. The 
short timescales behaviour is identical for both models and a simple argument allows 
one to rescale all simulations for different intermediate and final temperatures on the 
same leading curve which happens to be independent of the dynamics (identical for the 
symmetric and asymmetric KCIC models). 

We have also studied the FD relations for both models after a Kovacs quenching and 
a usual one. After the Kovacs quenching, the related FD relations verify surprisingly 
the FD theorem for both models in strong contrast with the results for the asymmetric 
KCIC model using the usual quenching procedure. No clear explanation has been found 
for this behaviour. It could be interesting to have a look at the Kovacs effect and the 
corresponding FD relations using the models introduced in [22] where analytical results 
may be expected due to the possible exact mean-field solution. 
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